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Abstract 


Linearized  supersonic-Flow  theory  is  employed 
to  evaluate  the  lift  and  drag  of  biplane  cellules  having 
the  "Busemann-biplane"  configuration.  The  lift  and  drag 
are  explicitly  expressed  as  functions  of  the  thickness  ratio 
and  the  angle  of  attack;  the  coefficients  involved  are 
universal  for  all  Busemann  biplane.  Interpretation  of 
the  results  for  various  Mach  numbers  is  afforded  by  a 
similarity  rule . 

Mast  of  the  results  are  presented  graphically. 

It  is  found  that  the  wave  drag  due  to  thickness  of  finite- 
span  Busemann  biplanes  is  small; 
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Introduction 


A  study  of  the  linearized,  equation  of  notion  of  super¬ 
sonic  flow,  with  the  velocity  potential  <p  as  dependent 
variable,  reveals  that  the  differential  equation  is  of  hyper¬ 
bolic  type*  As  is  well  known,  a  Cauchy’s  boundary  value 
problem  is  then  properly  set  if  initial  data  over  an  open 
surface  is  given.  This  is  exactly  the  case  with  a  supersonic 
wing.  Thus,  following  the  exhaustive  study  by  J.  Hadamard,  by- 
using  the  singular  fundamental  solution,  introducing  the  notion 
of  finite  part  of  an  improper  integral,  and  generalizing  the 
Green’s  formula  to  a  space  with  non-positive-definite  metric, 
a  solution  for  the  differential  system  is  obtained  which  bears 
exactly-  the  same  relation  to  the  distributed  boundary  value  as 
the  potential  to  its  corresponding  source  distribution  in 
Dirichlet’s  problem  in  the  elliptic  case.  Accordingly,  the  terms 
"supersonic  source"  and  "supersonic  distance"  were  introduced. 

Further  investigation  of  the  problem  shows  that  we  can 
always  introduce  a  fictitious  plane  separating  the  retrogressive 
characteristic  cone,  so  that  by  suitably  distributing  the 
boundary  value  on  the  two  surfaces  of  the  plane,  two  distinct 
solutions  can  be  obtained  for  the  two  separated  regions.  Thus 
the  introduction  of  a  diaphragm  and  the  exclusive  use  of  source 
distribution  for  the  determination  of  potential,  following  Eward, 
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can  be  justified.  The  detailed  proofs  of  this  "supersonic 
source-distribution  theory"  are  successively  developed  in 
Refs.  1,  2,  3o 

Thus,  the  study  of  supersonic  wings  resolves  itself 
essentially  into  the  same  study  of  the  potential  due  to  plane 
supersonic  source  distributions. 

In  the  linearized  theory,  the  Risemann  biplane  arrange¬ 
ment  becomes  the  one  shown  in  Fig.  1,  i.e.,  the  top  and  bottom 
surfaces  are  flat,  the  leading-edge  Mach  wave  of  either  wing 
intersects  the  other  wing  at  mid-chord,  and  the  airfoil  slopes 
are  related  by  the  formulas,  for  x  >  c/z 

Y/'x)  =-  ->i'  (x-  §) 

K  '(*)  -  -  Y,  '  ( x  -  J; ) 

The  typical  case  is  then  simply  that  of  two  isosceles  triangles 
pointing  at  each  other. 

In  this  investigation,  the  Busemann  relationship  between 
gap,  chord,  and  Mach  angles  shown  in  Fig.  l  will  be  assumed, 
but  it  will  not  be  necessary  to  specify  the  shape  of  the  profile 
in  deriving  some  general  results.  It  will  be  shown  that  the 
velocity  potential,  including  all  interaction  effects,  can  be 
calculated  by  means  of  integrations  involving  the  wing  surface 

slopes  oily.  The  general  results  will,  be  applied  to  the  - 

numerical  evaluation  of  the  wave  lift  and  drag  coefficients  of 
the  typical  Busemann  arrangement  having  triangular  wing  sections# 
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This  investigation,  conducted  at  Cornell  University,  was 
sponsored  by  Aeronautical  Research  Laboratory,  Wright  Air 
Development  Center,  Wright-Patterson  Air  Force  Base,  Ohio 
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Notation 


Coordinates 

Cartesian  x,  y ,  z,  ;  5,3,5 

ISich  j  j  2  j  5 

U  uniform  free  stream  -velocity  in  pc  -direction 
f  disturbance  velocity  potential 
pj  j3  pressure,  density  (  /*«  ,  J°fl  for  free  stream) 

M  free-stream  Mach  number 
/2>  -/ 

Cp  pressure  coefficient 
-n.  normal 

cr  local  slope  of  -wing  surface  in  flow  direction 
A  10cal  slope  of  diaphragm  in  flow  direction 

$  supersonic  source  intensity 
£  wedge  angle 
?<x  chord 
2.t>  span 

c  gap  =  V/» 

integration  area  over  surface  being  considered 
S'  integration  area  over  interacting  surface  of  other  wing 
Subscripts : 

T,  3  top  and  bottom  surface 
u.,  t  upper  and  lower  wing 
l,  zl  of  wing  and  of  diaphragm 
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PART  I.  THEORY 


A.  THEORY  OF  SUPERSONIC  BIPLANE 

'  (l)  Formula  for  Source  Distribution 

(2)  Calculation  of  Diaphragm  Distribution 

(3)  Solution  of  Integral  Equation 
(It)  Calculation  of  the  Potential 

B.  SIMILARITY  RULE  FOR  TIP  FLOW 

(1)  Conical  and  Pseudo  Conical  Flow 

(2)  Case  of  Busemann  Biplane 

(3)  Similarity  Rule 
(U)  Remarks 

PART  Ho  APPLICATION 
A.  POTENTIAL,  LIFT  AND  DRAG 

(1)  Potential 

(2)  Lift  and  Drag 
R>  COMPUTATIONS 

(1)  General  Form  of  Velocity  Potential 

(2)  Lift  and  Drag  in  Terms  of  Potentials 

(3)  Consequence  of  Similarity  Rules 

(li)  Lift,  Drag  Coeffid.  ent,  and  L/D  Ratio 
(5)  Computation  of  L  and  D  • 

Co  RESULTS  AND  CONCLUSION 

£ 
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PART  I.  THEORY 


A.  THEORY  OF  SUPERSONIC  BIPLANE 
Formulas  for  Source  Distributions 

The  equation  satisfied  by  the  disturbance  velocity 
potential  j>  in  the  linearized  theory  is 

where  subscripts  denote  partial  differentiation  "with  respect  to 
the  rectangular  Cartesian  coordinates  *  ,  y  ,  z  ®  Here  M 
denotes  the  fre e-stream  Mach  number,  and  the  coordinate  x  is 
taken  in  the  direction  of  the  undisturbed  stream.  It  has  been 
assumed  in  deriving  Eq.  (l)  $hat  9  4>«  ,  and  02,  are  small 
compared  to  the  stream  speed,  if  0  k  consistent  approximate 
formula  for  the  pressure  coefficient  is 
P  ~  Pa 

CP  =  S  pa  y*  =  ~2^y/v  (2) 

where  />,  ,  /j  are  the  pressure  and  density  of  the  undisturbed 
stream. 

An  elementary  solution  of  Eq0  (l)  is  the  so-called  super¬ 
sonic  source,  j>  (x,f,z)  = 

-/Z3'*'*  -  5  )x]  /z  ,  provided  that  the  value  zero 

is  taken  outside  of  the  Mach  cone  that  originates  at  the  point 
J  ,  3  ,  J  •  For  brevity,  we  shall  adopt  the  following 

notation:  _ ^ 

f*lz)  =  [  -  ^(y-D)1- p'z']  2 
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It  is  well  known  (Refs*  2,  3)  that  a  continuous  distribution 
of  these  singularities  over  a  surface  parallelto  the  flow  yields 
a  solution  satisfying  Eq*  tu  and  the  boundary  condition  2><%n  = 
on  the  surface*  Moreover,  Eward  (Ref*  li)  has  shown  how  a 
distribution  of  these  sources  over  a  fictitious  diaphragm  at 
a  wing  tip  can  be  used  to  account  for  the  interaction  of  upper 
and  lower  surfaces  of  a  monoplane  wing* 

We  shall  adopt  Evvard's  scheme  here  for  the  calculation 
of  tip  effects  for  both  upper  and  lower  wings,  placing  a  diaphragm 
at  each  wing  tip  and  introducing  the  conditions  that  these 
diaphragms  are  stream  Surfaces  of  the  flow*  The  potential  at 
points  on  the  top  ( r)  and  bottom  (£>)  surfaces  of  the  upper  (  u.) 
wings  is  given  by 

(3) 

-f  for/'tcjc'S 

s  s  (U) 

and  there  are  analogous  formulas  for  the  lower  f^wing*  The 
areas  of  integration  £  ,  on  the  wing  under  consideration, 
and  £'  on  the  other  wing,  are  shown  in  Fig,  2* 

Now  the  integrations  over  portions  of  £  and  S’  can  be 
simplified  immediately  by  use  of  monoplane  results.  First  of 
all,  it  is  dear  that,  in  all  areas  unaffected  by  biplane  inter¬ 
action,  the  wing-surface  boundary  condition  requires  that  £  —  U<r/n 
where  <r  is  the  slope  of  the  wing  profile  in  the  /<  direction* 
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Moreover,  Eward  has  shown,  that  fcr  monoplanes  -  and  therefore 
for  biplane  regions  unaffected  by  interning  interaction  -  the 
integration  over  the  diaphragm  can  be  replaced  by  another 
integration  over  part  of  the  wing.  For  any  point  forward  of 
mid-chord,  i.e.,  x  <  CL  ,  there  can  be  no  biplane  interaction, 
hence  it  is  convenient  to  write  the  relatively  simple  expressions 
for  these  points  before  going  on  to  treat  the  interacting  regions. 

%  <  cl  ;  no  biplane  interaction:  Here  monoplane  results  are 
applicable.  For  both  upper  and  lower  wings,  we  have 

K  (*,})  =  -  ^[<rrf4fo)aS  -  (5) 

**  *2. 

4>e  (*<])  =  ~  ¥J (6) 

X  >  cl  .  We  consider  now  a  point  on  the  upper  wing,  top 
surface »  If  the  point  lies  forward  of  the  Mach  line  from  the 
tip  mid-chord  (outside  of  area  a/  in  Fig.  3),  there  is  again  - 
no  biplane  interference  and  Eqs.  (5)  and  (6)  apply.  For  a  point 
in  a/  ,  however,  there  exists  an  effect  of  the  lower  wing, 

>> 

transmitted  through  the  interaction  regions  of  the  tip  diaphragm* 
We  can  write 

-%-(A«r'°>aS  (7) 

•Tr  u 

Here,  and  in  subsequent  formula,  s,  we  dendt  e  by  A  (  3  ,  ") )  the 
slopes  of  the  tip  diaphragms  of  upper  and  lower  wings.  Thus, 
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for  any  point  on  the  top  of  the  upper-wing  diaphragm,  ^  u  T 
is  equal  to  Tf  /\  u  ,  and  this  value  has  been  used  in 
eq.  (7)#  In  regions  unaffected  by  biplane  interaction  (e.g. 
for  ),  A  ( 3  ,  ?)  is  the  same  as  for  a  monoplane 

and  Eward’s  results  will  be  used  for  such  regions*  In  inter¬ 
acting  regions  A  is  still  unknown ,  of  course j  its  determination 
\ 

constitutes  the  main  problem  of  this  investigation.  We  shall 
postpone  this  to  the  next  section,  after  writing  an  analogous 
formula  for  points  on  the  bottom  surface  of  the  upper  wing® 

All  points  of  the  bottom  surface  of  the  upper  wing, 
for  which  X  >  are  affected  by  biplane  interaction® 

Let  Si  and  denote  the  areas  of  the  lower  wing  and  its 

diaphragm  that,  affect  the  point  .  The  wing-surface 

boundary  condition  is 

U&  <«'?)  -h^i4c[f*i7/'U>c'S  (8) 

Sr  Sjz 


This  is  an  explicit  formula  for  £  uQ  tx,#)  , 
involving  only  known  quantities.  It  may  be  noted  that  in  the 
region  S v'  >  t  lt  has  been  put  equal  to  V  A  l  / ® 

Moreover,  here  A  L  is  a  monoplane  value  unaffected  by  biplane 
interference,  and  is  therefore  known’  from  Eward’s  worko  We  now 
have 


J  r  +  ~>n 

-  ^f/Zn-fUcJoS'  -  %  f  A  LftOctS 

WADC  TR  52-276 


St 


S' 3 


(9) 


_q_ 


where  in  Si  is  known  from  Eq.  (8)  and  A  L  in  So 

is  known  from  monoplane  theory*  Again  the  calculation  of  the 
diaphragm  source  distribution,  in  So  ,  is  postponed 

to  the  next  section* 

For  the  lower  wing  there  are  formulas  exactly  analogous 
to  Eqs*  (7),  (8)  and  (9),  which  will  not  be  written  out  here* 


2.  Calculation  of  Diaphragm  Distributions 

The  conditions  that  insure  that  the  tip  diaphragms  will  be 
stream  surfaces  are  the  conditions  of  equal  slope  and  equal  pressure 
on  top  and  bottom.  Since,  as  Eward  has  pointed  out  (Ref*  5),  the 
diaphragms  of  a  rectangular  wing  tip  are  not  vortex  sheets,  equal 
pressures  imply  equal  values  of  f  ,  the  perturbation  velocity 
potential*  We  have,  then,  in  region  Sm  , 


d<pT  _  3  fa 
d  z  ~  aad 


(10) 


The  first  of  these  equations  leads  to 


7 fXu  C*'})  =  (11) 

=  -  Jus  '*'J>  Ji[ dr  i- 

C'  r.  J 


« 


r 
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The  second  Eq.  (10)  states  that,  in  Sit  9 

-  -  —j  Xu  fJ(0)<*S 

=  -J  t«3S(°)c*9-  7T f^LT -  wf*Lr(0°*£  (12) 

St  t  $n  St  $&< 

nhere  ^UQ  in  J1*  and  £&  is  given  by  Eqs*  (8)  and  (ll), 
respectively.  We  have  now  an  integral  equation  for  the  diaphragm 
slope  A  u.  t  for  points  7  ,  J  in  $ 


There  is  an  analogous  equation  for  XL  ,  which  will  not  be 
written  out* 

E<q.  (13)  is  to  be  satisfied  for  all  points  x  ,  y  on 
the  upper-wing  diaphragm,  -^or  some  areas,  there  is  no  biplane 
interaction,  i.e0,  jV  and  J1/  vanish,  so  that  the  second 
and  third  integrals  on  the  right  side  of  Eq.  (13)  disappear.  It 
is  clear  that  for  these  points  the  third  integral  vanishes  as  well, 
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since  Sr  and  S'n.  do  not  contain  any  points  5  *  ^  affected  by 
interactiono  Consequently,  for  non-interacting  points  ?  ,  g  , 
Eq*  (13)  reduces  to  E-ward* s  integral  equation  for  the  diaphragm 
slope  of  a  monoplane  (Ref*  ii)« 

2*  Solution  of  the  Integral  Equation 

Eq©  (13)  can  be  written  in  the  fora 


J =  F,  (x,  3) 


for  points  a  ,  v  in  ^  ,  where  2.  Ft  (x,  p)  denotes  the  entire 
right-hand  side  of  Eq.  (13),  and  involves  only  known  functions*  We 
now  introduce  the  new  coordinates,  ju.  f  \r  measured  al&ng  the  two 
families  of  Ifach  lines  on  the  wing  in  questions 

5  =  1  =  Fi 

HI  =  V  = 

|  I  £  (  **  i  A7  l  , 

?(*•>  =  [  'f*(7-Vx-/3'2']  "Z  (15) 

=  ((**-,  -*0  -  Mx7jX/^  j 

Our  integral  equation  now  takes  the  form 


I M,  - 


cX<r' 


lA7-  v 


—  F  (  ■u  I  j  ) 


for  points  Mi  $  ia  S jl  • 
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The  solution  can  now  be  found  by  means  of  the  following 
process : 


a 


We  now  multiply  both  sides  of  Eq*  (18)  by  1  / J  v-'~  v~t 

integrate  with  respect  to  ^  ,  and  exchange  order  of  integration 
in  a  manner  similar  to  that  just  employed.  The  result  is 
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■which  implies  (dropping  the  primes) 


This  solution  can  be  used  to  calculate  the  slopes  in 
regions  of  interaction.  This  completes  Eq.  (7)  for  f>ul  ,  and, 
by  use  of  Eq*  (11),  also  completes  Eq*  (9)  for  fu8  •  Eq*  (20) 
constitutes  a  generalization  of  Eward’s  egression  for  the  tip- 
diaphragm  slope,  to  Tiiich,  in  fact,  it  immediately  reduces  when 
-<4-  ,  v  lie  in  a  region  free  of  biplane  interaction. 


U*  Calculation  of  the  Potential 

Although  the  biplane  problem  is  now  completely  solved  in 
principle,  the  straightforward  calculation  of  ,  especially 
for  regions  of  biplane  interference,  ty  substitution  in  Eqs.  (7) 
and  (9),  is  extremely  tedious 0  Fortunately,  as  will  now  be  shown, 
it  is  possible  to  eliminate  entirely  the  integration  involving  X  ^ 
in  these  two  formulas . 

In  both  Eqs.  (7)  and  (9),  the  term  involving  At<.i» 


J  XlA['/o)*S 

Jr 


‘a  U  (*4,v)  d\S 
J  O',  -  \r 


(21) 


where  now  lie  in  region  Sz  • 
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We  return  to  Eq.  (13),  which  holds  for  points  in  S‘a$ 
and  write  it  in  the  fora 


*a 


'  rtu 

/J  M I  —  M 


Or  (u-jv,)  = 


Mn  ,//  , 

177  *  (m-'  v>) 


(22) 


where 


r^f  ,**  rs; 

=)  - Ij  -7jr=r-  -ft ] KTjJ-J  («y) 


(23) 


$'(*.,*)  =  {  f<rLTfu)c<S  -h  / Xtf4(Oc/S  { 

JS‘  Sa'  J 


(21) 


Actually,  (j> 1  is  the  potential  contributed  at  ^  ^  by  the 

lower  wingo 

The  solution  of  Eq»  (22)  can  be  written  down  immediately 
(Ref0  6)$  viz,. 


Cr(  <r, )  —  - 


M 


(25) 
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Since  Eq*  (22)  is  correct  only  for  points  m. ,  ,  v~,  in  fi¬ 

le,,  for  <  C~,  -we  must  restrict  ^  in  Eq*  (25)  as 
indicated* 

Now  for  the  points  outside  of  the  interaction  region,  i.e*, 
for  c  2  /'^  o',  ,  the  interaction  potential  '(u,\ 

is  zero*  Thus  Cr  C  <s->  ^  )  is  also  zero  for  <*<  * / 4  1/7  t 


We  can  now  consider  an  integral  involving  Cr  (^j  ^7)  j  i*e*, 

,  K-  .  n. 

I  k)  s  /  Gr{«yv,)  =  /  C?( m,  1/7) 

J  J M,  -M  J  J**,-**.  ' 


where  k:  £.  -u. , 


If  K.  <  ^  also,  &■(**- j  xr ^can  be  taken  from  Eq0  (25): 


alter  some  manipulation**  Recalling  the  meaning  of  Cr  ( Mj  \r, )  , 

(Eq*  (23)),  we  can  write  Eq.  (26)  as 

*  Ref.  3,  p0  36-UO* 

-16 
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j  K  /“/ 

J  if,  -f 


$  (“W>) 


Ik 

1  1 

n  [ 

k-u‘ 

M^c\ 


Ur, 


(f  M  ' 


7  4  f  *-**!  f  f  -J3li}?Lj.  ±[  4t<r  2>  ,  ,  ) 

3J^JZ¥^aj  J  Jr-.-v 


Since  the  only  restrictions  on  Eq.  (27)  are  lx  <  -'U.,  , 

and  k  <  v-;  ,  it  is  exactly  the  result  we  need  for 

Eq.  (21),  in  which  ix  “k  1/7  <  » 

We  are  now  prepared  to  write  complete  expressions  for  the 
potential  on  top  and  bottom  surfaces  of  the  upper  wing,  by 
substitution  in  Eqs.  (7)  and  (9)*  Let  Sj  be  the  portion  of 
Jr  for  which  -a  <,  u~,  ,  as  indicated  in  Fig*  3j  then 

=  -  wf<r^f'0jc,s  -  Tr,L(tr^-‘r^ro,c,s 

Sx  J  T- 


(28) 


-  -L-  f  yzt* («,*)<*$ 

JSZoi-Ta 

+  r~  /  -+-)*«' 

<?  w  J  J  -I f,  L  K  -M'  M,~  U'J 

-  ¥fjr"sro)dS  -  ¥  j  rr°> 

J*  Sr+Sx 

+  —  f  (<^B- <zT)f<fo)o/s  +  ^  f /^r°)  4z^('M^)dS 

2n  4^4 


C<j>' («•.*)  T^'[-L~  - —  j 0/u1 

zn  )  •  J  J  u,-f,  Lis-,-14'  u,- * -4’J 


-h 


MlCl/4  vt 

frx7^ 


(29) 
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Formulas  (28)  and  (29)  permit  the  calculation  of  the 
potential,  and  consequently  the  pressure  distribution,  on  the 
biplane.  It  is  seen  that,  whereas  we  have  succeeded  in  eliminating 
the  integrals  involving  \u  t  for  the  upper  wing,  we  are  left  with 
integrals  involving  A  L  ,  to  be  taken  over  certain  interaction-free 
areas.  In  fact,  if  interplane  interaction  of  a  higher  order  were 
encountered,  such  as  an  area  of  the  lower  wing  influenced  by 
interacting ' regions  of  the  upper  wing,  it  would  always  be  possible 
to  eliminate  the  \  integral  expressing  the  last  stage  of  tip 
interaction* 
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B.  SIMILARITY  RULE  FOR  SUPERSONIC  BIPIANE  THEORT 
(l)  Conical  and  Pseudo-coni cal  Flow  Field 

It  is  veil  known  that  the  flow  field  in  the  tip  region  of 
a  supersonic  flat  rectangular  wing  is  conical.  If  instead  of 
being  flat,  we  have,  denoting  by  cr  the  surface  slope  along  flow 
direction  0<  ,  that  c T  ex,  )  =  cr  (  50  only  -  i«>e,, 

along  each  chord  station  ;x  **  constant,  the  surface  slope  is 
spanwise  constant  —  then  we  have,  in  place  of  a  conical  flow  field, 
a  flow  field  which  is  "pseudo-conical,"  Thus,  denoting  by  <f>  the 
disturbance  velocity  potential,  the  conical  flow  field  will  be 


characterized  by 


f  =  **£(■%£) 


while  the  "pseudo-conical"  flow  field  will  be  characterized  by 


t  =  f 


where  =  J  M  z  —  /  ,  A7  ■  free  stream  Mach  number, 

While  the  conical  flow  field  finds  great application  in  the 
monoplane  supersonic  wing  theory,  it  will  be  shown  that  a  rectangu¬ 
lar  supersonic  biplane  flow  consists  in  addition,  in  the  inter¬ 
action  region,  of  a  "pseudo-conical"  flow  field.  Thus  the  success 
of  developing  a  biplane  wing  theory  depends  much  on  the  existence 
of  the  latter  type  flow0 
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Let  us  take  a  Btisenann  biplane,  and  follow  the  result  and 
notation  of  Part  I  where  the  potential  at  any  point  (?<■■/)  of 
upper  wing  top  and  bottom  surface  respectively  are  given  as 
followst 

tTN)  =  -  % -  £h*-„a-*-„r)r°>d* 

Si  Jj0 

-nL  C'yAre-xs. 

' w 

Ke1*'?)  =  ~  yffae/"0)0'*  ~  wj 

Jr  £*  + Sn 

^ 2n  I J  (J'L1B  (°)d$  ^  zn  f  '/*'*>&  PCu.'OdS 

Sz° 

r*  , -  (U) 


in  which 


CP  =  —  /T  —  —  oi 

U  J-  ° £.& 


S~u&  <SLr 


—  d  i~  <f 

—  -o/+  < T 


V  <  CL. 


X  p>  «_ 
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(5a) 


It  was  shown  in  Ref*  3  that 

a)  A  1  = 

TM  M  77  '  /W  7T 

.  Ujjt-y,  /  n  ,  . :  -•  ui  -  ^ 


IT  OiT  7  r/  (<£*> -tflrj  r 
L,  -  1 z 


I,  =-2 


z 

Asf  C 


'  *2.  /ft il+'/tfL-nxc'rI  * 

/  , ./  u,(u,r\r9)  -  m1cl/2_  ,  , - 

/  j.  —  ■  ~~  ydv^-a  _ 

'  ^  // A>/+  ^  V  *-  ^/,/f  cfrv,)^  +  A7  V  ^ 


/,= 


]  *4  ,*A  -  A7  V  — 


_/  U,!u,-r\f,)  -  M'-C'-/-l  (  ^  ■+ 

u.-v,  a-,  **>  -f~w'4 


J(ia,~v,)*-  '-/A'c  1 


v,*-1'/ r  , ,,  iv,-'/,  .  .  •  -/  2  ", 7 

f-  — /^tn  ,-—  =-=  =•  —  .-• — — - / 

2  ‘  ■J(u.  +  *Ji)  -M'c'-  a»  v  »-  J 


7T  ;  ,  ot,!--  l/,  '*'-  /■  l  <•/,  A,  '<  ' 

,  A/C  I  7  ) . . . 7 - - — rrr-- 

\  —  — —  I  1//7  L-  ‘f  +  '• 


•f 


,  _  /  tt.  ( tt ,  t-  V/)  -  M  'c  '/z. 


-Afe 


_/  ft< ,-  V)  )  -f  t  V 1. 


^ ■  J ( n,  r- - /h'-c'- 

J  ^  (  **  ,-i/J  -  ^  ^/2  -A  »0  7  /’•*/-i//j'Ly.  Aivf 


L  A  An  '-c  '-J 

j 

A 


*  {x-PP) 


Using  the  relation 

*'sZp'*+P?? 

P'Z)  =  { lX-5)*-^(p~0 


(5a)  can  be  put  in  form 

zr  * 


=  j3/  (*,  -ot)  (v,-^) 


k-  m **,(%.$) +  *ki*i>y] 


b) 


9^' 


f  7  =  Trfcrh  +(rLB-*lr)u] 


(5) 


(6a) 


7T  +  -/  _ _ _ 

1 3  ~  z  2A7  c  J  L4.  tr-t  _  A^xcV<? 


» 


77"  .  ,  /M  >  -t-  f^XC  '■/^  .  •  -/  ("l-t  </»  J  1 

/  *•  1  '  *  _ _  -f-  -  - - -~ 

F  «.J  «,  J 

,  -/ /m ^-i/7v;  -  /w  V1-  +  {^,  +  v,)  Jfu,~v,)l-+  *  v 
—  J — ■ — ■  —  -  - — 


J  /*4,  tv,)*-  Mxc^(  +  J  {O'.  -<S,)  V  '-z  ' 

2./W  C. 


■  Jy,  (u.-ia,)  -  Ai'cy^  + 


",  -r  j 


(  -ts,)  ^  -4  An  '< 
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(6a)  can  be  put  in  form 

2#/  U  f  ,  J  /l r  Me  rl  ,  vr  MC\] 


e) 


^  /a. 


^L*f£l  0^l4i 

-np  l  X  f 


■>  x  ,>  -I  P7  71  7 

x  +-T} 


-  -  *4'  -  ^  i  -  w(i-  $ ) 


3  / 

/  /o 


77^3  ' 

z/cnx 


if  (Fix’  [  .  -)  f>y  P'3  .  e-i  *  r  a  \  „  ■-<  P$/x 

~  ~Y  I  ^  -(/-f)*~  TZT/> ~ 

-  +  §^'^11  *  n./,. pj  ( ?a ) 


(7a)  can  be  put  in  form 


(TpttddS  =  ^  ir 

''j? 


a)  ,.j 

Jjro  *-a3  r*~3 

_  vt^e-Qr),  P  a"1 


ufats-Oir),  .  P*  ~  *k  p?  ■ 


1  *0-0) 
2  n  fl  /  (^C  ^T\  fa  ~  (Qua  ~  ^vr)  J f  J  [,i]  dj  J 


#-  fief.  3.  Appendix  5. 
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•  ^  +  $?<**'" W  ~  (/~  /I  */x) 

=  -  4sr^-'z^in  -  7~*(*  f2®  - ') 

*(•-$)  P?  *  1  *  (,~  %)  ' 

yjS(f^-J^)  *  ,(u/jjf,  -Kjagi  _$ 


so  (8a)  again  be  put  in  form 


-£fi 'r*a-W»S  -  ffM,&2)+f4r(4'V]  (8) 


Now  by  (6b) 


,  py .  vr  /•"/£?]* 

JruZ71  "/  L^-rr^>0^> 

erf*'  (?&?(£,<&)*(#) 
nM°  4  %.  $  ^Ja  J  %  -  m 

=  id,-  Tc)i  -■  4  Oj  ZfH) 

so  (6a)  takes  the  form 

-f„f (»:e-f„)r>^s  =^7/«'if (%'%)  +  siJei>%)\  .  , 

4.*J >  (6> 


U)  u, 


i.  (  b  ) 


Moreover, 


~  v^s,  [v>  ,  v,  wsJWzL_r_l _ 

=  ?£}£  [\ /j?  am  t  wii'  -g .  y±£ \!  ' 

nM  J 0  r  '  ^  y  P  77/^7  -  '  vt  ’  M<  *o-J\ 

'A,  !-Ma  0  1/7  _  fr-/3y  _*./,_#) 

u,  ,  +  P}/*  ^ 

~ 


WADC  'TR  52-276 


-23- 


“  v; 

Ij'fK  (£■  -  * ' = ZjrM,  (?■  9) '  sm  }>]  <?> 


(10a)  can  be  put  in  form 


f -  £XH  r#  s;  -  5j] 

Vr*4 


(10) 
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Substituting  (5)* (6),  (7),  (8),  (9),  and  (10)  into  (3)  and 
(U),  we  see  that 

i-  t-(^>  £  *•-*> kiS)  (ul 


4'  =  /j  3  ■  ^  •  denotes  regions  of  interaction,  which 
is  actually  a  "pseudo  conical"  flow  field. 


It  is  seen  from  above  that  in  eq,  (ll),  (12)  s  for  *  <  , 

a/x  term  does  not  occur,  so  flow  field  is  conical,  for  x  >  cc  9 
^/x  terms  occur,  so  flow  field  is  "pseudo  conical  ,n  However, 
for  both  regions,  the  same  similarity  rule  holds,  i.e,: 

"For  different  values  of  Mach  number,  differs  ty  a 

constant  factor  '/$  ,  provided  the  ratios  /3^/x  and  <l/^ 
are  maintained  constant," 

This  constitutes  a  "similarity  rule"  for  supersonic 
flow  at  different  Mach  numbers  o  It  leads  directly  to  the  determina¬ 
tion  of  a  Mach  number  correction  factor  for  converting  the'  biplane 
property  calculated  at  one  Mach  number  into  that  at  other  Mach 
numbers, 
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for 


iu  Remarks 

It  might  be  noted  that  the  above  Similarity  Buie  is  not 
limited  to  bi-plane  application  only.  The  same  holds  for  delta 
■wings  as  long  as  the  leading  edges  lie  either  outside  of  Mach  cone 
or  in  flow  direction,  and  O'  -  <r  (*)  .  However,  as  soon  as  the 

leading  edge  goes  inside  the  Mach  cone,  the  Siniiarity  Rule  ceases 
to  hold*  Indeed,  even  in  the  simple  case  of  conical  flow,  it  is 
easily  seen  that  Ts Chaplygin’s  transformation  will  give  for  each 
Mach  number  an  independent  boundary  condition  for  the  Laplace  equation, 

which  much  complicates  the  problem* 
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PART  n.  APPLICATION 


A.  POTENTIAL,  LIFT  AND  DRAG 
(1)  Potential 


With  the  Similarity  Rule  established,  the  stucly  of 
supersonic  biplane  is  greatly  simplified*  Indeed,  the  only  com¬ 
putation  needed  mill  be  for  the  case  of  M  —  Jz.  •  *or  all 
other  Mach  numbers,  the  required  conversion  relationship  uill  be 
furnished  by  this  Similarity  Rule* 


Using  the  general  expression  for  potential  derived  above. 


me  have 


-*s Is  ffhtyx 

t-,  =  -  %-Jsf‘.Tr'°J<x  -  if  -£rfe«~nr)/"‘W 


Sx 


-r„S  -  * 

JSoffz  ‘  '  ' 


^13  =  "  v~J +  Hf  (<rta  -  ^ 

/z 

-  rtjrv  jf‘*f  ?1'A:  -^1“" 
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■where  it  is  understood  that  4>/ ,  ^ ^'/sc.  311(1  &  come  in  only 
where  the  interaction  is  effected#  fy  carrying  out  the  integra¬ 
tion  as  given  in  Eef#  3#  we  have 


4  '=  -  ^Ti 

T M  M  77  2  t 


U-fdla-Glr) 


j  /  If  <3~*B  j  _  U  (Qua  -d hr )  j. 

^  —  wrr  '  A7  7T  z 


A*  _ 

T  —  p  [  f  (*r*)  -Mxc  %  _  j  Jcu, A^c  y+  7 

2'  <  /0  M,-U  *  Jo  **>  “J 

.  r?  v  .  /  ~  a» v*-  . 

=  r  ”E”  £  ~  ~ »  +v«<^--*r-  * 


'?c-  X*  «t<n ,+*!)- ,  "A  ZL 
2  ^  *  v 


- 


'/Vr 


2-  '  ^  Sfa+ijp-M' 


A »C  i-<  ^  Vi 

J  /ytA  y - — - 


^tj(*4(  +  Hxe' 


■  -/  ‘*AU,~''-> J 


Li,  ■f\T/ 


_  Mc  f  i,'-'  tAl±r±^Op^£yA  , . ;  -/  *-**  y 

-m)  -M-xy/4  ^({  +  j 

J  o  m,-<a  Ja  ■«,-« 


_  Mvc->-  A,  lf2 
-  ~  ~ 7  y 


<-OY-  <-0 


f/t1  ^7  F  .  •->  «!-  IS-,  .  .  -  -y  2r,  -J(t4,-\r,j’*'i-  a*1c  x  / 

— /  . . . —  :  .  -/-  — ,. ,  . -  . =—  / 

2  (  /A L,-t<yjA  mxc^  J  (**,  r  v->J  ^  -  a»  v  -»-  _y 


.  ty,''‘-/-f’'d'LA/A',y-  ^t)  J(u, -KiJ'  v-A’  V  ' 


^  J  (o,  i-  V,J  v-ssi'~cx  [  u,—*, 

--^W 2.  _  t-  -V  (M,  -  t^/J  A  A?  V  Va 


f  <f/J  -  tu'C  ‘ 


M  '-J(ui  ~<SiJ  '•y-'jV' 


?  c 


■+ 


(  &~LG  ~  fi'cr)  1%  J 

~  ^tr)  I2  J 


1/ 

// 


IL  y.  f  -1 

2.  ZA1C  Jut^_v'cv4 


- i 

2 


,  7j  ,  ,  _/  (h,-*;)  u,  t  M^Vz  ,  ,,  *  Vt  ] 

[  y  — /Un  - ;- . -  - ==---— 

**'  V fn 


**  >J~?  U ,  +  7  V  ■ 


.  1 -I  -t^lx'C~  t  (u,-h  vn  t  AlV-1- 

-  ...  .  _ _ _  —j  7-~-~-"-~— - — — 

V («.  «-/v,  -  <r,+ 


-t 


2  a/c 


Jn(u,-v,)  -  M^c/i  +  (*,  +n'c' 


“  —  'T*  4  +  . 


Tip  region  potential  of  a  symmetric  ning  for 
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X  <  a, 


(Appendix  5*  Ref.  3) 


ia  =  -  - 

rt  o  o 


z  >  *- 


■fu?-  *f 


(b)  *3 U  ffW  -Vr'J 


iff* ~-  r0i 


0  0 

i  f<r.  r*'Py 


-  f '*-*■> 

•  *. 


y-J 

v  v  "  F  v  (*&(* "  T 

(c)  /<.  =  "  iT ]  *J  j  <r,y(o)<fy  ~WJJ  +  P - *O/*,0)d1  +  V(<r;-rJ? 

=  tM- v-f 

-  ~^jP-p  -  B «  -  *(?£&&  **]*-*«; 

where  by  Appendix  5  *  •  3 

(i)  %  i-  y]^  =  W  *-  x^'ir  +  f  <*-*(*9) 

&)  J  +  2L]ctj  —  £  $  su^'iJ  -+  {3 y  y  ~  a, 

(3)  +  r°  -f2) 

-  rfc'Yf)  +  (*-*■) **^ '§^k  +  PI  ^  ~  t*- 
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o 


■when  #• 


2/<T/  _  g£gj> 

f<sO  X  <  0-  T  Ol.  ~  tj  2.  Zft 

<*>  !-*)-?*.(*-*)  =  %*- UJr  -  7a <«-°v 

<0  *-  -  }***  - 

-  Mi  U^a-ruTi /•"»<*$ 

jft  ...  XjJ-4 

J*a 


7/  rtf 

=  "i^r  JfJ  Ja 

=  ~znp(°u,a  '^tJ,  f  [  ££  Og-i)]s? 

~  ~ 2ng  ^8  -  °*tX  1/ 

u^.  _  </,_  ^  ,  t**r(?*-*  *1 

.  1^/2.  ~^T SZ  I  <*3  /  ,- . . . . - 


°i  A 


'rZ.fF-t'i 


~  “i/iyj  ~^r)zI,  ^ ^Tp  ( °Ze  ~°^t)J \  IL 

I, -Jj&Q  -*#(-**  -,)]* 

=  xsU?&-  +  ^''jy 

+  2*  izxZl'J  *‘py  -  -  2  Jfl-I*-?!'  - 

Ix  =  *£  +  P?  c*u~‘?j  -  1 *  M  ^  ’  77 

+  -frZjyf'')  ‘'zJfjtJrrt 

-  <LxU^'( -') 


♦Appendix  5*  Check  p«  H5jll£  Jtef.^g. 
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(2)  Lift  and  Drag 

With  the  velocity  potential  determined,  it  remains  a 
simple  matter  to  compute  the  pressure  coefficient,  lift  and 
drag#  By  the  linearized  theory,  the  pressure  coefficient  is 


given  as 


C,,  =i 


Af> 


2  u 

u 


_£  2$ 
V  I* 


a  Polf2- 

The  total  lift  and  drag  over  a  surface  are  obtained  from 


the  Cp  distribution  by  simple  integration 

With  linearized  theory,  the  lift  computation  can  further 
be  simplified  to  the  following  form,  which  es  especially  useful 
when  curved  surfaces  are  involved: 

z  -  =  ^(-d)ArA 

=■  -  Pa  zr  cty 

In  case  the  wing  has  an  isoceles  triangular  or  "diamond" 
profile,  since  the  constant  surface  slope  can  be  taken  out  of  the 
integral  sign,  the  determination  of  drag  simplifies  to  the 


following  form 


ft  V 


=  -  PoU  f  V  3  -  $  fj*)  ] 

For  a  multiplane  system,  with  projected  area  J1  in 
plan  form,  the  lift  and  drag  coefficients  of  the  system  are  obtained 

as  follows:  r  ^-(-0 

r  J 
_  &  h-  cr 

/ 

-32- 


C 
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B0  COMPUTATION  AND  RESULTS 


(l)  General  Forms  of  Velocity  Potential 

The  potential  at  a  point  tx,  )  or  (**■/  'f)  in  region  D  of 


the  upper  wing  is 


In  fact,  (1*1)  and  1.2)  are  general  forms  of  the 
potential  in  all  regions*  In  region  C  9  the  point  or 

considered  is  beyond  the  influence  of  the  rear  half  diaphragm, 
and  so  the  tem 

j_ 

Z7T 

vanishes  in  passing  from  region  _D  to  c  ,  and  remains  nil  on  the 
rest  of  the  wing  surface.  Further,  since  differs  from  zero 

only  for  J  on  the  wing  and  wr?  —  on  diaphragm,  the 
integral 
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goes  to  zero  as  <x,f)  goes  to  the  hyperbolic  boundary  separating 
regions  P  and  C,  •  In  fact,  in  all  the  regions  A#  B,  C  it 

<pu  /  (- 

contributes  nothing,  as  the  effective  area  for  ?  ^  in  *>!<,+  Ju 
disappears.  On  the  other  hand,  the  integral 


n 


degenerates  to 


-ijr  i-  Sjx 


-1 


fffo) 


77  Jfr 

in  regions  A,  B,  C,  In  the  regions  A  and  B#  i#ec,  those  areas 
forward  of  the  hyperbolic  a  lt  =  ,  the  effect  of  the  tip 

region  of  the  neighbouring  surface  cannot  be  felt  and  the  aja&ytic 
expressions  for  the  integrals  are  available. 

There  is  a  difference  in  the  surface  slope  c r  far  the  region 
3  <  cl  and  3  7  d-  9  but  the  formulas  (1,1)  and  (1,2)  remain  true. 

In  region  A^  the  tip  effect  of  the  supersonic  monoplane 
disappears,  the  flow  reduces  to  that  of  a  two-dimensional  Bisemann 
biplane. 


The  structure  of  the- present  theory  is  based  on  the  linearized 
theory  of  potential  flow.  As  a  consequence,  the  potential  will  be 
linear  in  the  angle  of  attack  d-  ,  and  the  thickness  ratio  ci,  in 
the  foim 


f»e  =  <r 


(1.3) 

(1.10 
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For  the  same  point  (x,£ )  on  the  loner  wing,  similarly,  we  have 

^T=  -  <&■*'*■  €*■  r 


-  tr  «  +  &■* 


(1.5) 


(1.6) 


(2)  Lift  and  Drag  in  Terms  of  Potentials 

In  the  Linearized  theory,  as  shown  before,  we  have 

+*■[<!%> 


(2.1) 


_  r  -,X,  f  r  f  J#/*'*- 

-Di/fu yw-l,  7  -  «j**Jx'f  *  */*<•  4,7  (2.2) 


where  the  (+)  sign  stands  for  top  surface  and  the  (-)  sign  for 
bottom  surface* 

With  the  distribution  of  in  the  Sisemann  biplane. 


°<,x  <  cu 

Q  <  x  <  2a 

-CX 

—si 

7~  c/3 

cT 

u-  <r 

CTt  T 

-x  +  r 

-  <j" 

(2.3) 


and  noting  that  <fl  f0 )  =  o  ,  we  have 

^  ur,  lz>  /pv  ~  •/  9*/T%a  i  <?'  J  #"*{2*0  * f 

J-UB  ,±-T  /fiir  —  -of./  f  J-.y 

^r,  Us/fiXf  =  o'*  f  <?%&«>  ±  otf-J (2 

■^uSy^-r ~/pU  —  ~  ^  •  Jcf‘0ie  tea)  ~  f  'J[2%e  (*)  - 

+  «sj[2tc)(a)-¥«3(Z*)+ 


(2  Ji) 
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She  total  lift  and  drag  of  the  biplane  becomes 


lAPU  =  *7 ~ 

JD/^U  =  «yf[ 

+  fXf  [  tV(Z^ 

q  =  f/fli f2 

L/gq  =  <*•  /  (  u  ?UT}fz*J  -  zr  t**  (*«)]<* 3 

=  o'*-/ / if  il*1^ -S- 

+  ^./ (§  *la<**}  -  tr 


(2.5) 


(2.6) 

(2.7) 

(2.8) 


3.  Consequence  of  Similarity  Rules 

a*  lift  and  Drag  over  the  Tip  Region 

Brom  the  Similarity  Rule  for  the  tip  flow  of  a  Busemann 


biplane 


rr 

Z/ 


± 
X  ' 


11 

x 


«,  O 


and  since  &  and  <f  can  be  separated,  we  have 


'iT 
U ' 


(J) 
-  / 


*r  //) 
^  *,■  ( 


±,£l- 

X  X  y 


where 


/  =  -T 

't  =  CXT,  u8,  L~Ty  L& . 


(3.1) 
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(3.2) 


Introducing  _T  =  ft  $  /2  a-  >  176  obtain 

4*  =  *  fi'S)  -iZ  <i.s>l<s 

=  f  L  ■* 


P^  =  *■(**)■  [i  -  */' *  v 


31 J 


(<rj 


?  -  ZJo  f  L*  a,*>  -to  (os)j^ 

;  =  iy 

o  1  tJT 


(«)  ;  ^  7  V  - 

r/J  (  ~.  SJJC  > 


(3*3) 


(3J0 


•where  l  and  p  are  universal  constants  for  all  ftisemann  biplanes. 

b*  Lift  and  Drag  for  Off -Tip  Region 

Beyond  the  tip  region,  the  flow  is  two  dimensional.  So 


(3.5) 


=  - 

77  «■ 

T5 

JL  d)'*',  G)  _  _ 
t;  /c(t3  1  G ~ 

7T  Q 

/* 

= 

I 

II 

•~\ 

rx 

S* 

Mb 

7T 

f 

a 

2«  -2" 

{*■ 

# 

and  we  have 


£■«!?"«)- IF  *17  i*)-° 


7TP,(«) 


Z  n <1 


(3.6) 


(3.7) 
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Let  L  o-H-  denote  the  lift  of  half  biplane  system 


(It*  5) 


l/  =■  ~ 
/ D  0/ 


1  *  zf4?+  Z  ~  4  ^ 


(lt*6) 


(5)  Computation  of  Z  and  £>  « 
As  shoTm  in  (3)  j 


^  =  n  J7u  '?'ur>(z«)-  S  *y 


2f 

5  ^  ±-  C  (  ?r  IL 

7TCL^  L  V 


>  r  7T  .  7T  .^7  7 

/ 7"  ^0  -  - 


(5.D 


(5.2) 


Setting  .  p*  —  '  j  ci  -  /  o  ,  i»e  have 


Z.  = 


/OO  7T 


/  0  0  7T  J  l  If 
0 


V  .  (<f)  rr  (<n  n 

77  91s  iz°)  ~  XT  2  Pua 


Let  ~j  crua  jt(o)  o/Z  =  -  /Lt  ■  <*  +  Af-  7 

l/  =•  &*'<*  +  S*‘S 

JZo 

=  Jv  ^  -.Zr  <T 


-*/,  A«>£|*W  = 

'/Jz  **  -Ta 


Vac 


Sr  Q-t-  Sj^  Sja  -Jit,  +  7 


Y  {z,n.«  -ff  (ril-e 

Sr  /  Tjl  Si  t  £tl 


Izhr 
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(5.3) 


(5.1t) 


(5£) 


where  an.  the  coefficients  of  o'  and  &  are  functions  of  (x;  )  • 

Through  (5.1) 5  the  general  forms  of  potential  (1*1)  and 
(l»2)  become  respectively 


Si.  +  f. 


For  regions  C  and  CB)Z  ,  we  have 

U  ft T  =  *  A  l^Aot  ~fzS*l  +  s'i  [  ° 

~T  =  dt  ±  \- JlAu+ 1  +  s-  -1  +Ji2r  7 

*  *  I'l  *kJ{  W]i  *l-T<-4JJ/d 

J  jr  *'  T 

For.  region  (A  )2  >  the  tip  effect  becomes  nil,  so 

U  ftr  ~  ^ '  'Ay  +  0 

Vts  “  '‘I -*****  ]  +f\*'-kJr  1 

i  -AiJJjjvU 


(5.8) 

(5.9) 


(5.10) 

(5.11) 
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For  region  (3Jt ,  which  differs  from  a  difference  in 

surface  slope  and  by  the  absence  of  interaction  between 
surfaces,  we  have 


C*ot  -B«]  +  & 

3L  fuB  =  ex'-  [-  +B«]  +  r[-*S+Qs] 


(5.12) 

(5.13) 


For  region  (A), 


iv-r  =  Qt-  +  <T-o 

=  ot-.j-As ]  +  *  f'A*l 


(5.U0 

(5.15) 


In  fact,  the  coefficients  A  d  ,  ,  A  <r  ,  -5^  *  Xf 

as  well  as  the  integrals 

JTvffj  7«'J  311(1 

can  all  be  analytically  evaluated.  Therefore  the  potentials 
in  region  A  and  &  can  be  put  in  puidy  analytic  form 

f  =  f”"*  t  f 

For  the  rest  of  the  integrals 

fs.  ,  niyi.  Jr*i,  /*■-*.*> 

Jr-tiz  Jzo+Jji 

numerical  computation  and  graphical  integration  becomes  necessary* 
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-V 


'g;c:  /|-s  n&z 


I2  = 


M,  -  ir,  ^  g.t  C4,+  i <7 _ 

^  J  (  *A  M^C 


~*7 _ -f^r'  ^  zd^^'l 


Ar£  r  ^ 


51  _  ^.-/  -rt*cK+{*t,+''i)J («.-*,)*-+  /wvc*- 

r t/O K-  v'x. ^[*4, - 1//  +  J/u.-\ 

_!_  ■<  (>~,u'(  *4r+'/0-*/)'cy_  Ai--i  *4,(0,-^))  +  m\'~/z 
**•  J(m<  +  'rOx-»i'c 1  M»  jsui-is,)  V/^  v'* 


/  "  Ol  ‘5 

V.-  •' —  ; 


(5.22) 


-  J ir,  («,-  io)  -  „\yz  -+  IT,  J f+"'cx‘ 

Note  that  the  limiting  value  of  <pu  'goes  to  zero  along 
the  hyperbolic  boundary  off  the  mng  and  thus  9^1s 

continuous  across  the  hyperbola.  On  the  other  branch  of  the 
hyperbola  on  the  "wing,  '  satisfies  the  relation  (Bef •  3  Appendix  1;) 


ft  ■ '  -  -  v~fe(y -f‘>  -  -  *  -J  +  ^  <- 


•which,  in  fact,  is  the  expression  for  the  potential  $/for  region  4 
and  43  *  These  remarks  serve  as  a  check  to  the  numerical  computa¬ 
tions* 

Since  in  region  C  ,  X>  and  over  the  diaphragm  the  expression 
for  PJ  is  complicated,  we  carry  out  the  computations  numerically 
for  points  in  these  regions*  As  ^/only  enters  into  the  integra¬ 
tion  over  the  regions  extended  by  the  inverse  Mach  cone,  we  need 
only  compute  for  points  bounded  by  the  hyperbola  mv  —  , 

the  trailing  edge  of  wing  and  the  inverse  Ifach  line  on  the  dia¬ 
phragm  extending  from  trailing  edge  tip*  Hence,  using  /2  =  /  , 
cl  —  C  —  / o  ,  the  range  of  computation  will  be  as  follows* 
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«  IT  = 

VA--  2.5,  3  6  ,  7,  9.  '/,  /j 

*4  •  X  7,  <9,  9y  /o,  /!,  12 J  /3>  /J",  /  7,  )S  ,  2 !  ,2  1.  23,  £4 

Our  computation  starts  with  a  chosen  (say  =  it),  and  pro¬ 
ceeds  -with  the  formulas  (5.21),  (5*22)  in  computing  1,  and  I2  for 
a  range  of  xt* 

Figo  1;  &  5  present  the  results  for  I,  and  1 2.  •  Fig. 

6,  7  give  X-  (u,  *■)  and  (u,*). 

For  the  sake  of  convenience,  we  use  the  principal  branch 
in  the  identity 

<Scn-X  *4.  —  ( u  j-  J  u  '~-i ) 

?or  determining  cX  */  ^ 


(5.2)  Computation  of  ty  c.  (for  region  D  and  C  ) 

5  57 

=  °ItV  +  (^a-VC7)lx' 

—  -7o4  •  ^  +  Jj  ■  S’ 

Where  =  //_  21J 

7<r  *  2/ 

By  Appendix  3.  Ref.  3: 

7  /  77  ,  /  ^  7  -  ^  1 

3/  =•  —  7-  - ~ - 


(3.23) 

(5*24) 

(5.25) 

(5  .26) 


2-M  C  J 

**  '  j  (Lf,  -  ^ 1  y-  A7  V  u  A7  V  *- 

—  /U^T1  yl>  ^  /  M!  +■  'S'jJ  fA4‘~  vOZ  A  ^  v  c  L 

„  c  -^'^XA  -/  Vf]  fut-i/i)  y- 


(5.27) 


M'C  ^  [  C4,-  V,  +  Jb,-  </,)'-+ 


(5.28 
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9  £2,’ 

The  limiting  value  of  along  the  hyperbola  mit-jq  is 
zero  off  -wing  and  constant  —U  <£~7on  the  -wing*  so  it  i»  continuous 
across  the  curve  ^  a~=  o  •  is,  ho-wever,  discontinuous 

across  the  midchord*  All  these  have  been  used  as  a  check  for  the 
numerical  result* 

The  numerical  data  and  the  computation  procedure  are 
similar  to  those  for  ^  *  In 


J /-t  y 


—  fcf!  v 


the  principal  branch  of  x  is  used.  In  fact,  much  of  the 
computed  result  for  cfjcsxi  be  taken  over,  and  so  the  labor  of 
calculation  much  reduced* 

Fig*  8,  9  present  the  result  of  ,v)  and  7 f  \j)  * 


(5*3)  Graphical  Integration  of 


ffj*}  and  f  fft) 


•4  V, 


Since  potentials  in  regions  c  and  D  enter  into  the  expressions 
for  lift  and  drag  only  as  ( (z  cl)  ,  we  are  concerned,  therefore, 
only  with  those  points  (m.,  ,  m  J  along  the  trailing  edge  pr  -  'z 

For  numerical  confutation,  we  choose  the  following  four  points 

=  4  /  3,  2  <3  4,  3),  (/r.2S-4jH)j  ( /5-.  z8<^ '3J)t 

(  S0Jl,  Jo/1  ) 

For  the  last  point  at  tip,  the  limit  of  the  integral  converges  to 
(Ref*  3)* 

^Joi  ( 'Of?,  /of)-  V  -  -rrj&  (  ,ajit  >oJi).  <T 
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Writing  the  integral  in  the  form 


A*  '  !  ’  I  k 

j  1  A  j 


-  7  •  [  V.  /  ,  r  ,  /  ^  '  f  ,  ,  I 

Jso/r,  U--w 


//^  i  J  M  i  ~  i/~i 


~  Z  *  V/v,  -  ^ 


f  Ji  ( "  ’s'/,)  OUi' 

Lo/  u>)j v, -«> 

^  (5.2 9) 

L  u.-»wT^  W-w 

'K/'  -f-  V,  /  | 

^3“^  ~f,’  (^-e,  K)  J 


€  =  o.  /a,  we  get 


t'  A/-y  /  -r°/^ 


-  r^*' J*'  {u'.irjtfu1  ^  , 


v'- *w 


(5.30) 


The  error  involved  will  be  of  the  order  6  2  =  a.  <?  o  J  # 

After  the  point  iyr  y  was  chosen,  we  simply  have  to 


form  the  integral 


v/y  \  ^ ) 

Cn ,-w) 


(A-'-  Of,  <fj 


over  She  interval  f  ,  \rt  -  o ./ o)  on  the  diaphragm  off  wing, 

with  {  <a’  j  v-,  )  and  ^,)  known  from  (6*1).  The 

graphical  integration  is  performed  by  using  a  planimetero  The 


last  correction  term 


is  approximated  by  . 


The  result  is  plotted  in  Fig0  10. 
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(5.U)  Graphical  Integration  of 


'  and 

tLctv 


It  is  seen  that  the  second  integral  is  contained  as  a  part  of  the 
first  integral,  and  can  easily  be  separat  ed  from  the  first  at  the 
last  step  of  integration*  The  values^'  and have  been  computed 
in  (6.2)  9 

The  integrals  are  evaluated  graphically  for  points  (  -*<.  , ,  ) : 

(I0j2,  10j2.) t  (I5.28ii,13),  (17.28^,11),  (19.28U,9)  in  region  D, 
and  (21.28i*,7),  (22.28U,6),  (23.28i*,5),  (2l*.28U,li)  in  region  C 

for  which  the  second  integral  disappears  and  the  first  reduces  to 

f  f  Ja<,  J  d* 


f  ^ M-  f  J**' ( 

Jn  J  J  J  v, 


In  each  of  the  stripwise  integration  with  respect  to  v-  and  4c  , 
the  same  technique  as  explained  in  (6*3)  is  adopted  for  the 
neighborhood  of  the  singular  point*  Results  are  plotted  in  Fig.  11,12 


(615)  Computation  of  —  j  (7^  a  ^  ( o)  ot  ^ 

fz 


As  shown  in  Part  1^  A(l),  ( i )  j 
For  region  (B)g  and  C 

Jr  0 

—  -  j-  [  ^ A  -t-  (&-x  ]■  w 

(  (5.31) 

_  JL  f  -l  a  +  (&-x  +  {??)  -( x  ~ajJ  ■  f 

n  L  n 


For  region  D 
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(5.32) 


where 

A 

a 


~fln  3  -fill  f 

- 

—  [ ’-<**£''  yy]  P7 1  *  ' 't  *  t  (x-  zfyJ 

=  *•  -ZJ'J 


C  =  A 


o 


j  r~r~  ‘  ^  r-  sT~~  ^ _ 

and  (4,  ^  =  g*  v-  d  >  ^cy#  —  cJ  -  <f  are  noted* 

Along  f  =  o 

-  f£*r"»*t-fl-*Tr*+=K-]s 

Qr  putting  cT  —  °  and  changing  the  sign  of  d*  in 

we  get  -  f  C£r  (4  ( o) d$ 

Jfi  1 

The  points  ^  )  coaputed  are*  ^  -  20j  #  -°»  1*615, 
lio)44i,  7.272,  10.101,  11.515#  12.929,  and  lli.3U3»  The  results  are 
shown  in  Fig.  13. 


(5.6)  Computations  of  =r 

Pro.  Part  Itt(l),  (ii),  nth  <rJz,-cru'P=  2dt  C,  cr^/-  %‘K  i  x-f, 

we  have: 


for  region  B^  ,  C 

i  /mb -fur)r° •>  1, 

J/"o  / 
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for  region  D 


^Z°  r  7  /  T  ^ 


=.  •  o<  /-  ^Z'  <r 

t  fi 


2  -  -I 

Z  Z  L  *-/ 


where 


I,  =  tttz'p^p- 1-0} -2* 

4  =  * ^  2 -'  +f}[a-t-<Z 

+**[■ i a'Jif,  - -U2 Jff,} -*Jft f- 'f*7#  -J*'-a-f?) 

The  same  set  of  points  on  the  trailing  edge  as  in  (6.5)  is 
used.  The  result  is  shown  in  Fig.  llu 

(5*7)  Computation  of  Potential  in  Region  and  Bj) 

The  potentials  outside  of  tip  interaction  region  (C  and  D) 
are 

tT  - '  ff’^T-Sfajcts' 

Lb  (5.35) 

te  =  -JU.r^S'  -  £,//<”> far -5-  -  4  ' 

JI  s^b 

*”  </.<  -  Zfr*T ro*i ;  <5*36) 

(see  Part  IA.(8))» 


(5.36) 


jo  _f/4 


i)  Ibr  i.e.,  /, '  =  J/ =  o  ,  so 

9  —  44 _j-  r O  W  /  .  ri 

^f/3  ~  tt  ~  Ti  'di-o) 

=  £  t-d  +  i) 

*  ~[ f 


(5.37) 
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If 


[  f /4 folds].  ✓  -  [f  [/(o)  <f '].  cT  (5.38) 

ha  Jx»  hi. 


Now,  with  ]■  -  / 


ro^y  =  Tip,  *  zJjfJTj  -  2  *  tf^'J 


if  /*'v*S  =  ~  11  * 

^  77  * 

4-  I*  -h  ■£  X  -<£uo  /— 2 
^  X 


(5.39) 


(5.1*0) 


So  ln  U  ^  if  ?'•  lot  +  ~4‘,  'S>-'T  It'-u7,uji)  ( S'  4lJ 


v  1*1 


xj  t  r  =  ft*™  ^  U  fjp  -  -  iff*^ 

C  -  ff/f  *r  W  K’  =  -j£Ht  - 


-  -//’  re> *!  v  ?«?  =  '  ry ”"f 


(SJii) 


^  etc,  are  in  analytic  fonns0  i  o 


ii)  For  i.e„,  o.<  v<  20.  ,  f  y 


2/  =  71  l/=  0 
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■where,  with 


J  f*(o)  d  S  -  ^  JJJ  x- y  -  Jx/UZf 

J.  J  /4/o)ef$  =  -JjjZj  +  £ y  c^r'j 


7T  .  /  w  £ 

“  ^  ~x 


4^'J  =  v  a.  -*■  -  *  *  Z^'/^Z- 

/  fMVoip  —  77  /'•jr-^J 
. r' 

So  J  and  ^/J^are  again  analytic* 


So  /  UQ  and  are  again  analytic* 

The  limiting  values  of  these  integrals  on  the  boundary 
separating  A  and  B  are 

//^  X  =  77*  ^  [nflJ  ffeJo/S'=  rra^ 

Jra  j/->x  JrXa 

/J/oJotf  =.  -  7T  fy-a.) 

y->iJsfh  7-**  Jrr' - 

which  check  the  computed  results* 

In  the  computation,  the  integrals  such  as  (5*37),  (5*38) 
etc.  are  first  evaluated,  then  'f  M  r  '  /4-J  etc*  The  computations 

are  made  for  points  (y,  3O  as  follows: 

for  ff  2a)  :  *  = 

$  —  J 3  00  /  <2&Jz.-2o,  ^  y  Jz  —7  o  j  2? o . 
for  cp  (0l)  ;  o(  -  /6 

J  -  o,  !.(,/?,  7.2  72. 

The  point  ( 20  y  J^ao  )  is  the  intersection  of  hyperbola  aa  50 
and  the  trailing  edge  X  =  2.  a-  =  2  o  *  Results  are  plotted 
in  Fig.  15,  16. 


(5*8)  Computation  of  Potential  For  and  k^ 
For  A-^t  =*** 


(5.1*5) 


Z  j>  s  -  -  (0/+  5)f^Uo)<t$  -  -  77  x  lol+  s) 
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(5.1*6) 


For  A2* 


ZL  <S) 
ZJ 


-  qc  .  f  p'OJ  °*r  =  -rrjrc* 

Jfr 

=  -f'"*JIus“f  -f 

<r/°  vi,v  T;-' 

-  [  f  ^ro)  +  f  I > 

L  Jfx"j  Jfr'  J 


(5.U7) 


-(  2  a~y)  -TT  &  -  -n  2  <T 


(5.1*8) 


vr  Zon  tr  /)  <<2)  .  t  -n  ~0*) 

U  tfue.  (d)  ~  rr  ^ ur  ^  ~  u  r 


—  -  77  4- 


7T  /, 

-  <?^r  «j  =■  o 


=  -2a7r 


(5.U9) 


(5.50) 


(5.9)  Resume  of  Results  of  Confutation 

The  results  of  computation  for  potentials  in  region  D  and 
C  are  combined  according  to  (5.6),  (5.7),  (5.8),  (5.9),  (5.10) 
(5.11)  etc.  to  give  the  final  form  of  potentials,  i.e. 

Tt’  ( X )  —  ft-  fx)  •  0/  9^'  ^  ^ 

where  ^ ’  =.  /ty  tt,  ^  ^  ;  of  —  <3- ,  2  <a_ 

The  resultant  etc.  are  platted  in  *ig.  17. 

Finally,  the  functions 

and  .j~  jcTj  , 
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C.  RESULTS  AND  CONCLUSION 


From  the  result  s  of  computation,  it  is  obtained  that 

Z  =^2.3  27 

Jp  —  o  ■  -?  //  3 

Thus  the  lift  and  drag  coefficients  of  a  Bxsenann  biplane  with 
chord  length  2a>  span  2 bt  wedge  angle  <T  at  Mach  number  M  , 
and  angle  of  attack  are  given  by  (li.U),  (I»*5)  and  (U*6)  as 
follows : 


CP  =  7*5?  '  *I~3)  * 1 '  **  *'J 

=  ^  (2) 

2 /■  =  ±  _ L _ 

/p  <*  /  /  0-323  fS.,  (3) 

•7/3^  -33  -^6  1  V  J 


where  ^  =  /w2-y  ,  /??  = 

G*V 

Figure  (19)  gives  t  a  plot  for-£"<-  (which  is  also  /<x) 

against  Mach  number  M  with  different  aspect  ratioso  Figure  (20) 

„  (c ri)  / 

gives  the  same  plot  for  Lo  /  S z  .  Thus,  for  any  given  aspect 
ratio  of  the  biplane  at  a  certain  Mach  number.  CL  and  Cp  can 


Figure  (21)  gives  a  plot  of  L/D  against  p/R  for  different 


cx  ratios* 
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The  limiting  case  of  the  present  theory  occurs  when 
i#e.,  {3  <42.  —2.  ,  in  which  case  the  two  tip  Mach  cones  on  same 

wing  meet#  These  limiting  cases  are  either  given  in  plots  ty  the 
point  (g>  ,  or  indicated  with  incorporated  tables o  It  may  further 
be  mentioned  that  these  limiting  cases  lead  to  a  unique  result 
for  L/D  ratio  as  follows t 

(  z./\  —  1  _ L _ 

1  &  / -h  °/7i  ( r/ol)^ 

From  the  computed  result,  it  is  easily  seen  that  a  Busemann 
biplane  of  finite  span  highly  approximates  a  flat  plate,  as  the 
wave  drag  due  to  wing  thickness  is  inherently  low# 
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Value  of^/vant  % 

Limiting  Values  hiMitmtj  j 


>fH£ 


co%OfWoi 


/2 
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0,58/7 

7 

2.  236 

A/64 

2 

/.4/4 

2,327 

4 

/.//& 

4.654 

d 

A  ©<90 

6,980 

e 

/■031 

9,308 

/o 

AOZO 

II.  64 

/2 

1,014 

/J.  96 

-  Q/v  - 

_4 _ 
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-  v<x  — 

Im*~i 
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